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Abstract 

We propose a collective Hamiltonian which incorporates interactions capable to 
generate rotations in nuclei with simultaneous presence of octupole and quadrupole 
deformations. It is demonstrated that the model formalism could be applied to 
reproduce the staggering effects observed in nuclear octupole bands. On this basis 
we propose that the interactions involved would provide a relevant handle in the 
study of collective phenomena in nuclei and other quantum mechanical systems with 
reflection asymmetry correlations. 

The properties of nuclear systems with octupole deformations [0J are of current interest 
due to the increasing number of evidences for the presence of octupole instability in 
different regions of nuclear table 0. Various parametrizations of the octupole degrees of 
freedom have opened a useful tool for understanding the role of the reflection asymmetry 
correlations and for analysis of the collective properties of such kind of systems 0, |], 
As an important step in this direction it is necessary to elucidate the question: which are 
the collective nuclear interactions that correspond to the different octupole shapes and 
how do they determine the structure of the respective energy spectra? The physically 
meaningful answer could be obtained by taking into account the simultaneous presence 
of other collective degrees of freedom, such as the quadrupole ones. 

In the present work we address the above problem by examining the interactions that 
generate collective rotations in a system with octupole deformations. Based on the octa- 
hedron point symmetry parametrization of the octupole shape {§], we propose a general 
collective Hamiltonian which incorporates the interactions responsible for the rotations 
associated with the different octupole deformations. It will be shown that after taking into 
account the quadrupole degrees of freedom and the appropriate higher order quadrupole- 
octupole interaction the model formalism would be able to reproduce schematically some 
interesting effects of the fine rotational structure of nuclear octupole bands. The study is 
strongly motivated by the need of theoretical explanation of the recently observed stag- 
gering patterns in octupole bands of light actinides as well as by the possibility to gain 
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an insight into the fine structure of negative parity rotational bands based on octupole 
vibrations. 

Our model formalism is based on the understanding that the collective properties of 
a physical system in which octupole correlations take place should be influenced by the 
following most general octupole field V 3 = X)L-3 a 3% Y£ , (i n the intrinsic, body-fixed 
frame) which can be written in the form [Q: 



^3 = eoA 2 + £ei(^i« +E £ 2(W) , (1) 



i=l i=l 



where the quantities 



i , 1 105 

A = -^(Y 32 -Y 3 ^) = -^—xyz, (2) 

= no^lg^-^-^), (3) 

^1(2) = -jV5(y 3 3-n-3) + ^V3(y 3 i-^3-i) 

= \J^x(x 2 - - y 2 - - z 2 ) , (4) 
^i(3) = -i-Vbtyzs + Ya^-i^VziYu + Y^) 



F 2 (l) = ±(Y 32 + Y 3 „ 2 ) = ^^z(x 2 -y 2 ), (6) 

F 2 (2) = 1^/3(^33 ->3-3) + JV5(r31-K3-l) 



= 4\/^^ 2 -- 2 )' (7) 
r d V Iott 

f 2 (3) = -^^3(r 33 + y 3 - 3 ) + ^^5(^31 + ^3-1) 

' ' (8) 



r 
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(with r 2 = x 2 +y 2 + 2; 2 ) belong to the irreducible representations (irreps) of the octahedron 
group (O). A 2 is one-dimensional, while F\ and F 2 are three-dimensional irreps. The seven 
real parameters e and e r (i) (r — 1, 2; i — 1, 2, 3) determine the amplitudes of the octupole 
deformation. Their relation to the a 3 ™ is given in [|J. 

Our proposition is that the general collective Hamiltonian which incorporates the 
shape characteristics of the octupole field (fj) can be constructed on the basis of the above 
octahedron irreps. For this purpose we introduce operator forms of the quantities A 2 , 
Fi(i) and F 2 (i) (i = 1, 2, 3) in which the cubic terms of the Cartesian variables x, y and 
z in Eqs ©-© are replaced by appropriately symmetrized combinations of cubic terms 



2 



of the respective angular momentum operators I x , I y , I z (with I 2 = I 2 + I 2 + I 2 ). The 
following Hamiltonian is then obtained: 



with 





2 3 

H oct = H A2 + H Fr{i) ' 
7*=1 i = 1 


(9) 


H An 


1 . . ... ... 

= a2~ \[I<rIn + LJt)It + I z \I<rIn + i«ir)l i 
z ^ L v x y 1 y ^ / z ] z \ y 1 y x / 1 ' 


(10) 
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H F 1 (2) 


= ^/i 2 (54 3 -34/ 2 ) , 


(12) 


Hfi(3) 


= ^/i 3 (54 3 -3y 2 ) , 


(13) 


Hf 2 (i) 




(14) 


Hf 2 {2) 


= f22{iJ 2 -il-iJl-i%) , 


(15) 


H F 2 (3) 


= f 23 (i y i 2 z +i 2 J y + i 3 y -i y i 2 ) 


(16) 



The Hamiltonian parameters a 2 and f r i (r = 1, 2; i = 1,2,3) are formally related 
to the parameters in ([]]) as follows 02 = 60^/105/(4^), = ei(i) y^7/ (47r), ^2* = 

e 2 (i)^/l05/(167r), i = 1,2,3. 

During the procedure described above, the r 3 factors appearing in the denominators 
of Eqs are replaced by J 3 factors. In the final result, Eqs (P^|) -(|l6|), we normalize 

with respect to J 3 , i.e. we multiply the results by J 3 , an operation which is equivalent to 
the transition to a unit sphere, a natural thing to do since we are interested in surface 
shapes. 

We remark that the terms of the Hamiltonian obtained (as a function of the angular 
momentum operators I x , I y , I z ) correspond to the same octupole shapes which appear 
in Eqs (0)-© and belong to the same irreps of the octahedron group. In other words, 
through the above procedure we determine the octahedron point symmetry properties of 
the system in angular momentum space. 

Our analysis shows that the operator Hfi<i)i Eq. (|Tl~D, which corresponds to Y 30 (with 
axial deformation) is the only one octupole operator possessing diagonal matrix elements 
in the states with collective angular momentum /. Below it will be shown that it is of 
major importance for determining the fine structure of collective bands with octupole 
correlations. Actually, it is well known that the Y 20 (axial) deformation is the leading 
mode in the systems with reflection asymmetric shape (See for review 0). 

Further, it is known that the use of the pure octupole field ([[]) is not sufficient to 
incorporate the collective shape properties of the system. More specifically a unique 
parametrization of the pure octupole field in an intrinsic frame has not been obtained 
yet in a consistent way 0. In this respect the consideration of octupole degrees of 
freedom together with the quadrupole deformations is important. A general treatment of 
a combined quadrupole-octupole field is proposed in the framework of a general collective 
model for coupled multipole surface modes 0, |8j . 



Based on the above consideration we suggest that the most general collective Hamil- 
tonian of a system with octupole correlations should contain also the standard (axial) 
quadrupole rotation part 

H rot = AP + A'f 2 , (17) 



where A and A' are the inertial parameters. In addition the following higher order diagonal 
quadrupole-octupole interaction term (corresponding to the product Y" 2 o • ^30) could be 
introduced: 

H qoc = f qoc ^(l5i 5 z - UI 3 J 2 + 3lJ A ) . (18) 

This operator is normalized with respect to the multiplication factor I 3 . (More precisely 
we use the product I 3 Y 2 o ■ Y30) so as to keep all non-quadrupole Hamiltonian terms of 
the same order.) 

Then the Hamiltonian of the system can be written as 

H = H bh + H rot + H oct + H goc . (19) 

Here 

H bh = H + fj z , (20) 
is a pure phenomenological part introduced to reproduce the bandhead energy in the form 

E bh = E + f k K , (21) 

were E and f k are free parameters. The /^-dependence of E b h, which can be reason- 
ably referred to the intrinsic motion, provides the correct value of the bandhead angular 
momentum projection K in the variation procedure described below. 

We remark that the Hamiltonian (|19|) is not a rotational invariant in general. It 
does not commute with the total angular momentum operators and any state with given 
angular momentum is energy split with respect to the quantum number K. Therefore, 
the physical relevance of this Hamiltonian depends on the possibility to determine in an 
unique way the angular momentum projection. The basic assumption of our consideration 
is that K is not frozen within the states of the collective rotational band. We suggest that 
for any given angular momentum it should be determined so as to minimize the respective 
collective energy. The resulting energy spectrum represents the yrast sequence of energy 
levels for our model Hamiltonian. We remark that similar procedure is used in Refs. 



TO I in reference to the AI = 2 staggering effect in superdeformed nuclei. 



As a first step in testing our Hamiltonian we consider its diagonal part 

H d = Hbh + H rot + H oct + H qoc ■ (22) 

were the operator H^ ct = Hf^i) represents the diagonal part of the pure octupole Hamil- 
tonian H oct , Eq. (|). 

The following diagonal matrix element is then obtained: 

E K (I) = E + f k K + AI(I + l) + A'K 2 + f n (^K 3 -^KI(I+l) 

+ f qoc ^(l5K 5 -UK 3 I(I+l) + 3KI 2 (I + l) 2 ) . (23) 
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Following the above assumption for the third angular momentum projection, we de- 
termine the yrast sequence E(I) after minimizing Eq. (|23|) as a function of integer K in 
the range —I < K < I. The obtained energy spectrum depends on six model parameters: 
Eq essentially responsible for the bandhead energy; /& which provides minimal energy for 
K = K bh = I bh ; A and A' are the quadrupole inertial parameters which should gener- 
ally correspond to the known quadrupole shapes (axes ratios) of nuclei; fn and f qoc are 
the parameters of the diagonal octupole (PI) and quadrupole-octupole (p~8|) interactions 
respectively. We consider the latter two parameters as free parameters. 

We applied several exemplary sets of the above parameters and obtained the corre- 
sponding schematic energy spectra. One of them is given in Table 1. It is seen that 
the "yrast" values of the quantum number K gradually increase with the increase of the 
angular momentum /. We remark that they correspond to the local minima of Eq. fl23|) 
as a function of K. This is illustrated on Fig. 1. We see that these minimums are well 
determined and their depth increases with the increase of the angular momentum. Such a 
behavior of the spectrum corresponds to a wobbling motion and could also be interpreted 
as a multiband-crossing phenomenon. The obtained yrast sequence can be considered 
as the envelope of the curves with different values of the quantum number K as it is 
illustrated in Fig. 2. 

In addition we see that the K- values of the odd and the even sequence of levels are 
grouped by couples which imply the presence of odd-even staggering effect. Indeed, the 
presence of such an effect is demonstrated in Fig. 3 (a)-(e), where the quantity 

Stg(I) = 6AE(I) - AAE(I - 1) - AAE(I + 1) + AE(I + 2) + AE(I - 2) , (24) 

with AE(I) = E(I+ 1) —E(I), is plotted as a function of angular momentum I for several 
different sets of model parameters. (The quantity Stg(I) is the discrete approximation of 
the fourth derivative of the function AE(I), i.e. the fifth derivative of the energy E(I). 
Its physical relevance has been discussed extensively in Refs || [TT|.) 

Fig. 3(a) illustrates a long AI = 1 staggering pattern with several irregularities, which 
looks similar to the "beats" observed in the octupole bands of some light actinides such as 
220 Ra, 224 Ra and 226 Ra M. Also, it is rather similar to the staggering patterns observed in 



rotational spectra of diatomic molecules [12|. In Fig. 3(b) the increased values of fn and 



f qoc provide a wide angular momentum region (up to / ~ 40) with a regular staggering 
pattern. The further increase of f qoc results in a staggering pattern with different ampli- 
tudes, shown in Fig. 3(c). These two figures resemble the staggering behavior of some 
rotational (negative parity) bands based on octupole vibrations [O . The further increase 



of fn and f qoc leads to a staggering pattern with many "beats", as shown in Fig 3(d). 
Notice that in Fig. 3(d) the first three "beats" are completed by I « 40, while in Fig. 
3(a) the first three "beats" are completed by / « 70. An example with almost constant 
staggering amplitude is shown in Fig. 3(e). It resembles the form of the odd-even stag- 
gering predicted in the SU(3) limit of various algebraic models (see Ref. |§ for details 
and relevant references). It also resembles the odd-even staggering seen in some octupole 
bands of light actinides, such as 222 Rn ||. 

Now we can discuss the general Hamiltonian structure fll~9|) including the various non- 
diagonal terms (10), flT^) — Q16Q ■ Here, the major problem is the circumstance that K 



is generally not a good quantum number. However we are able to provide our analysis 
for small values of the respective parameters which conserve K "asymptotically" good. 
This requirement assumes a weak i^-bandmixing interaction which guarantees that for 
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any explicit energy minimum appearing in the diagonal case the corresponding perturbed 
Hamiltonian eigenvalue will be uniquely determined. Thus we are able to obtain respective 
K-mixed yrast energy sequence. Our numerical analysis of the Hamiltonian eigenvector 
systems shows that the parameters of the non-diagonal terms should be by a order smaller 
in value than the parameter / n . In addition, we established that the following couples of 
non-diagonal terms give the same contribution in the energy spectrum: H^ 2 and Hp 2 ny, 

H Fl (2) and Hp^y Hf 2 (2) and Hp^y 

In Fig. 3(f) a staggering pattern with a presence of .fT-bandmixing is illustrated. In 
fact we added the following three non-diagonal terms Hp^y ^Fa(i) and Hp 2 i 2 ) to the 
already considered diagonal Hamiltonian (^), with the parameters of the latter being 
kept the same as in Fig. 3(b) (and in Table 1). We see that the mixing leads to a 
decrease in the staggering amplitude with the increase of angular momentum so that the 
staggering pattern is reduced completely in the higher spin region. This pattern resembles 
the experimental situation in 218 Rn and 228 Th |J (odd-even staggering with amplitude 
decreasing clS cL function of J). 

So, the staggering patterns illustrated so far (Fig. 3) cover almost all known AI = 1 
staggering patterns in nuclei and molecules. The amplitudes obtained for the examined 
sets of parameters vary up to 300 keV. Some reasonable theoretical patterns with Stg(I) ~ 
500keV can be easily obtained. On this basis we suppose that the model parameters can 
be adjusted appropriately so as to reproduce the staggering effects in nuclear octupole 
bands as well as in some rotational negative parity bands built on octupole vibrations. 
Also, an application of the present formalism to the spectra of diatomic molecules could 
be reasonable. 

Here the following comments on the structure of the collective interactions used and 
the related symmetries would be relevant: 

1) The equal contribution of the three couples of non-diagonal terms (mentioned above) 
indicates that only four octupole Hamiltonian terms are enough to determine the energy 
spectrum. This result reflects the circumstance that in the intrinsic frame three octupole 
degrees of freedom, from the seven ones, are related to the orientation angles. For example 
we could suggest that the following terms [applied in Fig. 3(f)] give an independent 
contribution in the total Hamiltonian: Hp^iy Hp^), #f 2 (i) and Hp 2 (2)- We remark 
that our analysis (related to the collective rotations of the system) gives a natural way to 
determine the four collective octupole interaction terms. 

2) From symmetry point of view we remark that the diagonal term Hfui) which 
corresponds to Y 30 possesses an axial symmetry while the non-diagonal terms Hp 1 ^), 
Hp 2 {i) and Hf 2 {2) (of previous item 1)) are constructed by using the combinations (Y31 — 
F3-1) with C 2v symmetry, (Y" 32 + Y3-2) with T d symmetry and (Y33 - F 3 _ 3 ) with D 3h 
symmetry. So, our analysis shows that the axial symmetric term should play the major 
role in the structure of the collective rotational Hamiltonian while the non-axial parts 
could be considered as small i^-band-mixing interactions. From microscopic point of 
view, a detailed analysis of the above spherical harmonic combinations and the respective 
symmetries has been provided on the basis of the one particle spectra of the octupole- 
coupled two-level model 



3) The observed influence of the non-diagonal Hamiltonian terms on the fine structure 
of our "schematic" spectra suggests an important physical conclusion: the non-diagonal 
K- mixing interactions suppress the staggering pattern. In such a way we find that the 
axial symmetric term Hpux) is the only one pure octupole degree of freedom which provides 
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"beat" staggering behavior of the quantity Q2~3| ) (See Fig. 3(e)). (The quadrupole-octupole 
term H qoc gives an additional contribution and provides wider angular momentum regions 
with regular staggering.) So, our analysis suggests that the AI = 1 staggering effect 
observed in systems with octupole deformations could be referred to the dominant role of 
the axial symmetric "pear-like" shape. 

In addition, it is important to remark that the fine (staggering) behavior of our 
schematic energy spectra reflects the structure of the interactions considered through 
the K- sequences generated in the above minimization procedure. Thus our analysis sug- 
gests that in the high angular momentum region some high K band structures should be 
involved. From microscopic point of view the values K — 0, 1, 2, 3 have been included in 
calculations, showing that in the beginning of the rare earth region the values K = 0, 1 
are important for the lowest 3~ state, while in the middle of the region the values K = 1, 



2 are important and in the far end of the region the values K = 2, 3 are important [15 



The same authors deal with nuclei with A > 222 in Ref . |15j . One of the authors of Refs. 



IB] , [16] in Ref. |17| finds that the restriction to K < 3 is not justifiable for large energies. 
This is in agreement with our findings of Table 1. 

In conclusion, we remark that the collective interactions considered in this work sug- 
gest the presence of various fine rotational band structures in quantum mechanical systems 
with collective octupole correlations. In particular, they provide various forms of stagger- 
ing patterns which appear as the results of a delicate interplay between the terms of pure 
octupole field and the terms of high order quadrupole-octupole interaction. The analysis 
carried out outlines the dominant role of the axial symmetric "pear-like" shape for the 
presence of a AI = 1 staggering effect. The obtained multi K- band crossing structures 
could be referred to a wobbling collective motion of the system. We propose that the in- 
teractions involved would provide a relevant handle in the study of collective phenomena 
in nuclei and other quantum mechanical systems with complex shape correlations. 
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Table 1: The "yrast" energy levels, E(I) (in KeV), and the respective K- values obtained 
by Eq. @ for the parameter set E = 500keV, f k = -7.5keV, A = 12keV, A' = 6.6keV, 
f u = 0.56keV, f qoc = 0.085keV. 



I 


E(I) 


K 


I 


E(I) 


K 


/ 


E(I) 


K 


1 


522.772 


1 


13 


2335.81 


5 


25 


5453.12 


11 


2 


568.327 


1 


14 


2576.57 


6 


26 


5694.49 


12 


3 


637.095 


1 


15 


2827.57 


6 


27 


5935.5 


12 


4 


728.71 


1 


16 


3082.36 


7 


28 


6157.5 


13 


5 


840.857 


2 


17 


3344.94 


7 


29 


6378.29 


13 


6 


971.155 


2 


18 


3608.18 


8 


30 


6575.37 


14 


7 


1123.22 


2 


19 


3877.05 


8 


31 


6770.62 


14 


8 


1288.09 


3 


20 


4143.16 


9 


32 


6937.23 


15 


9 


1472.71 


3 


21 


4413.03 


9 


33 


7101.62 


15 


10 


1668.56 


4 


22 


4676.45 


10 


34 


7232.21 


16 


11 


1880.56 


4 


23 


4942.01 


10 


35 


7360.44 


16 


12 


2101.68 


5 


24 


5197.18 


11 


36 


7449.45 
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Figure Captions 



Figure 1. The diagonal energy matrix element Ek{I) (in MeV), Eq. (p3|), is plotted as 
a function of K for / = 1,2,. ..,10, for the parameter set E Q = 500keV, fk = — 7.5keV, 
A = 12keV, A' = 6.6keV, f n = 0.56keV, f qoc = 0.085keV. 

Figure 2. The diagonal energy matrix element Ek{I) (in MeV), Eq. (E31), is plotted as 
a function of / for K = 10, 11, 12, 13, for the parameter set of Figure 1. 

Figure 3. AI = 1 staggering patterns [Eq. (pl|) 1 obtained: (a) - (e) by the diagonal 
Hamiltonian (^) for several different sets of model parameters; (f) by adding three non- 
diagonal terms H Fl ( 2 ) [Eq. (0)], ^f 2 (i) [Eq- (0)] and Hf 2 {2) [Eq. (p~5|)1 to the diagonal 
Hamiltonian (p2|). 
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